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Abstract: We study Kapitza-Dirac diffraction of a Bose-Einstein con- 
densate from a standing light wave for a square pulse with variable pulse 
length but constant pulse area. We find that for sufficiently weak pulses, the 
usual analytical short-pulse prediction for the Raman-Nath regime contin- 
ues to hold for longer times, albeit with a reduction of the apparent modu- 
lation depth of the standing wave. We quantitatively relate this effect to the 
Fourier width of the pulse, and draw analogies to the Rabi dynamics of a 
coupled two-state system. Our findings, combined with numerical model- 
ing for stronger pulses, are of practical interest for the calibration of optical 
lattices in ultracold atomic systems. 

OCIS codes: 020.1335 Atom optics; 020.1475 Bose-Einstein condensates; 020.1670 Coherent 
optical effects; 050.1940 Diffraction 



References and links 

1. C. Adams, M. Sigel, and J. Mlynek, "Atom Optics," Phys. Rep. 240, 143 (1994). 

2. P. Meystre, Atom Optics (Springer- Verlag, New York, 2001). 

3. P. L. Kapitza and P. A. M. Dirac, "The reflection of electrons from standing light waves." Proc. Cambridge Phil. 
Soc. 29, 297 (1933). 

4. D. L. Freimund and H. Batelaan, "Bragg Scattering of Free Electrons Using the Kapitza-Dirac Effect," 
Phys. Rev. Lett. 89, 283602 (2002). 

5. P. L. Gould, G. A. Ruff, and D. E. Pritchard, "Diffraction of atoms by light: The near-resonant Kapitza-Dirac 
effect," Phys. Rev. Lett. 56, 827 (1986). 

6. P. Berman, Atom Interferometry (Academic Press, San Diego, 1997). 

7. M. Inguscio, S. Stringari, and C. Wieman, (eds.) Bose-Einstein Condensation in Atomic Gases, Proceedings of 
the International School of Physics "Enrico Fermi" Course CXL (IOS Press Amsterdam, 1999). 

8. Yu. B. Ovchinnikov, J. H. Miiller, M. R. Doery, E. J. D. Vredenbregt, K. Helmerson, S. L. Rolston, and 
W. D. Phillips, "Diffraction of a Released Bose-Einstein Condensate by a Pulsed Standing Light Wave," Phys. 
Rev. Lett. 83, 284 (1999). 

9. J. Hecker Denschlag, J. E. Simsarian, H. Haffner, C. McKenzie, A. Browaeys, D. Cho, K. Helmerson, S. L. Rol- 
ston, and W. D. Phillips, "A Bose-Einstein condensate in an optical lattice," J. Phys. B: At. Mol. Opt. Phys. 35, 
3095 (2002). 

10. J. H. Huckans, I. B. Spielman, B. Laburthe Tolra, W. D. Phillips, and J. V. Porto, "Quantum and Classical 
Dynamics of a BEC in a Large-Period Optical Lattice," http://arxiv.org/abs/0901.1386 (2009). 

11. J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper-Kurn, D. E. Pritchard, and W. Ketterle, "Bragg Spec- 
troscopy of a Bose-Einstein Condensate ," Phys. Rev. Lett. 82, 4569 (1999). 

12. S. Gupta, K. Dieckmann, Z. Hadzibabic, and D. E. Pritchard, "Contrast Interferometry using Bose-Einstein 
Condensates to Measure h/m and a," Phys. Rev. Lett. 89, 140401 (2002). 

13. G. K. Campbell, A. E. Leanhardt, J. Mun, M. Boyd, E. W. Streed, W. Ketterle, and D. E. Pritchard, "Photon 
Recoil Momentum in Dispersive Media," Phys. Rev. Lett. 94, 170403 (2005). 

14. S. Inouye, A. P. Chikkatur, D. M. Stamper-Kurn, J. Stenger, D. E. Pritchard, and W. Ketterle, "Superradiant 
Rayleigh Scattering from a Bose-Einstein Condensate," Science 285, 571 (1999). 

15. D. Schneble, Y. Torii, M. Boyd, E. W. Streed, D. E. Pritchard, and W. Ketterle, "The Onset of Matter- Wave 
Amplification in a Superradiant Bose-Einstein Condensate," Science 300, 457 (2003). 



16. P. J. Martin, B. G. Oldaker, A. H. Miklich, and D. E. Pritchard, "Bragg Scattering of Atoms from a Standing 
Light Wave," Phys. Rev. Lett. 60, 515 (1988). 

17. S. Gupta, A. E. Leanhardt, A. D. Cronin, and D. E. Pritchard, "Coherent Manipulation of Atoms with Standing 
Light Waves," C.R. Acad. Sci. Paris 2, 497 (2001). 

18. W. R. Klein and B. D. Cook, "Unified Approach to Ultrasonic Light Diffraction," IEEE Trans. Sonics and Ultra- 
sonics 14, 123 (1967). 

19. O. Morsch and M. Oberthaler, "Dynamics of Bose-Einstein condensates in optical lattices," Rev. Mod. Phys. 78, 
179 (2006). 

20. I. Bloch, J. Dalibard, and W. Zwerger, "Many-body physics with ultracold gases," Rev. Mod. Phys. 80, 885 
(2008). 

21. R. J. Cook and A. F. Bernhardt, "Deflection of atoms by a resonant standing electromagnetic wave," Phys. Rev. 
A 18, 2533 (1978). 

22. H. Batelaan, "Colloquium: Illuminating the Kapitza-Dirac effect with electron matter optics," Rev. Mod. Phys. 
79, 929 (2007). 

23. D. Pertot, D. Greif, S. Albert, B. Gadway, and D. Schneble, "Versatile transporter apparatus for experiments with 
optically trapped Bose-Einstein condensates," http://arxiv.org/abs/0907.1323 (2009). 

24. M. Cristiani, O. Morsch, J. H. Muller, D. Ciampini, and E. Arimondo, "Experimental properties of Bose-Einstein 
condensates in one-dimensional optical lattices: Bloch oscillations, Landau-Zener tunneling, and mean-field ef- 
fects," Phys. Rev. A 65(6), 063612 (2002). 



1. Introduction 

The diffraction of matter- waves from a standing light wave is a fundamental concept in atom 
optics (HE). Originally predicted by Kapitza and Dirac for electrons more than 75 years ago 
(and recently also observed for these (4)), it was first demonstrated in the 1980s with an atomic 
beam [5 ], and has since become a standard tool in atom interferometry for the coherent mixing 
of momentum modes (6] [2). The advent of Bose-Einstein condensates in the 1990s (7) has 
made it possible to directly observe the dynamics of matter-wave diffraction in time-of-flight 
images GO |9] [TO). Moreover, the diffraction of condensate atoms from standing light waves 
finds applications in high-resolution spectroscopy ifTTTl and metrology (T2l [T3l, and plays a 
fundamental role in superradiance |[T4l[T5ll . 

Two diffraction scenarios can generally be distinguished. In the resonant Bragg case, the 
atoms oscillate between two resonantly coupled plane-wave momentum states depending on 
the strength and duration of the interaction with the light field lfT6l . If the Bragg condition is 
not met, atoms can nevertheless be diffracted into a number of off-resonant momentum states, 
provided that the interaction is sufficiently short and strong Q. In accordance with a common 
convention in atom optics ifTTl we refer to this case as Kapitza-Dirac diffraction, but we extend 
it to include times beyond the Raman-Nath regime. In this context, it is interesting to note a 
similar discussion for the diffraction of light from sound waves in acousto-optic devices |[T8l . 

In the Raman-Nath regime, i.e. when atomic motion during the interaction with the light field 
can be neglected, the populations of the diffracted states exclusively depend on the product of 
the strength Vb and the duration T of the interaction, i.e. the area of the applied pulse 0. 
Outside the Raman-Nath regime, the diffraction dynamics exhibits collapses and revivals for 
constant interaction strength (HI [9] [T0]|. In this paper, we now specifically analyze the case of 
a pulsed interaction of variable duration but constant pulse area. This allows for a study of the 
breakdown of the Raman-Nath prediction and, in particular, for a quantification of deviations 
when the system is close to, but not deep into, the Raman-Nath regime. 

The findings of our study are of direct interest for the elimination of systematic errors in 
calibration measurements for experiments with ultracold atoms in optical lattices [19, 20 ] when 
the lattice depth is determined via Kapitza-Dirac diffraction. 

This paper is organized as follows: Section [2] reviews general aspects of the matter- wave 
diffraction from a one-dimensional optical lattice, while Section [3] briefly describes our experi- 
mental system. In Section|4]we analyze diffraction patterns for weak pulses, and give an analyt- 



ical modification to the Raman-Nath diffraction formula, which is motivated by a comparison 
of the diffraction dynamics with that of a coupled Rabi system, as well as by considering the 
spectral properties of the pulse. Section [5] discusses the more complicated dynamics of strong 
pulses, and Section[6]deals with calibration methods for the depth of optical lattices, comparing 
numerical simulations with single- shot diffraction patterns not restricted to the Raman-Nath 
regime. 

2. Raman-Nath regime 

We first briefly review general aspects of the diffraction of a condensate from a standing light 
wave with wavenumber k = 2n/X that is switched on for a duration t. The standing wave 
gives rise to a sinusoidal optical potential Vbcos 2 fe lfT9l . The evolution of the condensate in 
the standing wave can then be modeled (2TJ [22 ] (neglecting mean-field interactions) as that of 
a matter wave y subject to the Hamiltonian 

H = -(h 2 /2m)d 2 +V cos 2 kz (1) 

where m is the atomic mass. Expanding the condensate wave function in the basis of plane 
waves populated by diffraction from a standing wave as \j/(t) = Y.n c n{t) e l2nkz (where n = 
0, ±1,±2, . . . and c n (t = 0) = S n ^) and introducing the dimensionless parameters 

a = (EP/h) t (2) 
P = (Vo/»)t, (3) 

(where = (nhk) 2 /2m denotes the ^-photon recoil energy, with E^ = E r ), transforms the 
time-dependent Schrodinger equation into a set of coupled differential equations 

dc n an 2 P , 

1-7- = C n + -7- (c n _i+2c w + C n+ i) (4) 

at t 4t 

for the amplitudes c n (f) of the diffracted orders n. For a given lattice depth, the highest momen- 
tum order (±2nfik) capable of being populated is given by the cutoff 

n = ^/Pj^c (5) 

for which the potential energy is fully converted into kinetic energy (cf. [10 ]). Note that in these 
equations, a is the pulse duration % in units of the 2-photon recoil time ffl = h/E?\ and j8 
measures the area of the pulse. 

The dynamics of the condensate in the standing wave depends on the ratio between an 2 and 
j3, i.e. between the kinetic energies acquired during diffraction and the depth of the potential. 
The Raman-Nath approximation consists of neglecting the an 2 terms in Eqs. which is 
justified if T is much shorter than the harmonic oscillation period in a potential well [8], such 
that zcOhnj^ 1, with C0h o = [VoEr 2 ^] 1 / 2 /h, or equivalently j3a <C 1. In this case, the solution 
of Eqs. <Q is c n (t) = (-/)^~^ /2t /„(/^/2t), such that the population P n = \c n \ 2 of the nth 
diffracted order after application of the pulse is given by 




where the J n are Bessel functions of the first kind. 



3. Experimental procedure 

In the experiments described in this paper, we subject an optically trapped 87 Rb Bose-Einstein 
condensate to a vertically-oriented, pulsed standing light wave at 1064 nm, for which the 2- 

(2) 

photon recoil time T r v 1 « 20 /is. A description of our apparatus for the production of conden- 
sates is given in ref. l23l . In brief, in a nearly isotropic crossed-beam optical dipole trap we 
produce condensates typically containing 5 x 10 5 atoms in the \F = l^mp = —1) hyperfine 
ground state without a discernible thermal fraction. The standing light wave, with a Gaussian 
l/e 2 radius of 130 /im at the position of the condensate, is derived from a single-frequency 
ytterbium fiber laser with a line width of 70 kHz and can be switched off within 1.5 /is using 
an acousto-optic modulator. Immediately after application of the pulse, the atoms are released 
from the optical trap and are imaged after 15 ms time-of- flight using near-resonant absorption 
imaging on the cycling transition. 

A series of images taken for pulses of variable duration t but with constant pulse area j3 = 4.5 
is shown in Fig.[T] Based on a naive invocation of Eq. ([6]) for the Raman-Nath regime one would 
expect the same diffraction pattern in all the images. Instead, the number of diffracted orders is 
seen to generally decrease with the duration of the pulse, with an oscillatory decay of the first 
diffracted orders for long pulse durations. 
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Fig. 1. Time-of-flight (TOF) absorption images of a condensate diffracted from a 1064 
nm standing- wave optical pulses of constant area /3 = Voz/h = 4. 5, and durations ranging 
from t = 4 ^s (a = t/t^ « 0.2) to 360 jUs (a « 18). For durations % > 75 jj,s (for 
which [P I a] 1 / 1 < 1.1), the apparent modulation depth of the standing wave undergoes an 
oscillatory decay consistent with the form (j3/2)sinc(a/2) (see text). 

The populations of the Ohk and ±2fik momentum orders in the series of Fig. [T] are shown 
in Figs. [2] (A) and (B), together with corresponding data for a weaker standing wave with 
P = 1.5, for which diffracted orders 2nhk with \n\ > 1 are not populated over the full range of 
pulse durations. For the stronger pulse with j3 =4.5, this is only the case for pulse durations 
exceeding 75 /is due to contributions of the ±4hk orders at shorter durations, as shown in 
inset (C). 

4. Weak-pulse dynamics 

For sufficiently weak pulses with a cutoff n ~ 1 for which only the lowest orders ±2hk are 
populated, Eqs. ^ can be reduced to three coupled equations; this approximation is valid for 
shallow potentials Vq < 4E r . The solutions cq and c±\ to Eqs. |4]), which can be obtained in a 



straightforward way, then lead to 



p±i = wr^ sin2 y r 2 " j (7) 

and Po = 1 — 2,P±i . This can be cast into the standard form of a Rabi oscillation 

1 /2 

with generalized Rabi frequency £1 = [x 2 + A 2 ] , resonant coupling # = Vb/ a/2^ and detun- 

ing A = Er /h = COr between the atomic momentum states Ohk and ±2hk. Another way to 
interpret this result is the following: each scattering event that changes the atomic momentum 
from to ±2hk by energy conservation leads to a frequency mismatch of the scattered photon 
by A with respect to the standing wave, requiring a corresponding Fourier width of the pulse for 
stimulated scattering to be able to occur. Consequently, the drop-off of populations in higher 
orders with increasing a can be seen as resulting from the decrease in the Fourier width of the 
pulse, which motivates a more general modification of Eq. ^ to 



Vi8 2 /2 + a 2 



P„=4^sinc|J, (9) 

where the sine function arises from the Fourier transform of the square pulse. Indeed, the pop- 
ulations Po,±l predicted by Eq. ^ agree with those of Eq. ^ up to 0(a 2 /3 2 ), independent 
of the pulse duration a or pulse area j3. For [/3/a] 1 / 2 < 1, the agreement of Eq. ^ with the 
experimentally observed patterns is excellent, as can be seen in Fig. [2] (A) and (B). 

It is interesting to note the connection of Eq. ^ to the resonant case of nth order Bragg 
diffraction from a moving standing wave (with a relative detuning 8 CO = n 2 A between the two 
beams). In this case, the Fourier transform of the applied potential Vo cos 2 (kz + S CO t/2) with 
duration T would lead to a frequency distribution oc sinc[(ft> — Sen) t/2] (i.e. identical to that of 
the standing wave case, but shifted by 8 CO) making available the resonant frequency component, 
independent of the pulse duration. 

The form of Eq. ^ shows that although the populations P n in a spectrum may follow a 
Bessel distribution as predicted by the Raman-Nath Eq. ([6]), a fit will not necessarily return 
the correct modulation depth. The apparent pulse area in Eq. ^ is given by (J3/2) sine (a/2), 
which means that for a pulse of duration T, the apparent modulation depth of the standing wave 
extracted from a diffraction spectrum is reduced from Vb to 

V ,eff = V sinc(T/2T r (2) ). (10) 

The apparent pulse area returned by a fit of Eq. ^ to the diffraction patterns is shown in Fig. 
|2|(D ) for both the j3 =4.5 and /3 = 1.5 data sets. For the latter, the fit results agree with Eq. 



P . a 



(10) over the full range, consistent with the fact that orders \n\ > 1 are not populated. The first- 
order revivals can thus be interpreted as arising from a sampling of the Fourier transform of the 
pulse. For j3 = 4.5, the observed behavior still shows qualitative agreement with revivals from 
sampling of the sinc-shaped frequency distribution with increasing pulse duration T, but clearly 
exhibits deviations when higher-order momentum states are present. 

5. Strong-pulse dynamics 

For more intense pulses with n > 1, the simple generalization of Eq. ^ no longer holds, as is 
evident by the large discrepancies in Fig.[2]for T < 75 /is in the j3 = 4.5 data set. 




Fig. 2. Suppression and revival of atomic diffraction from a constant-area standing wave 
light pulse with increasing pulse length. (A and B) Relative populations of the central con- 
densate Ohk and orders ±2hk vs. normalized pulse duration a = ijir ~ t/20 /is. Filled 
blue dots and open red squares refer to pulses of area /3 = 4.5 (cf.Fig. 1) and 1.5, respec- 
tively, and the solid lines plot the function J% [(/3/2) sine (a/2)\. The dashed black lines cor- 
respond to numerical fits for j8 = 4.5. (C) Decay of the ±4hk orders in the /3 = 4.5 data set, 
with the numerical simulation (dashed line) as well as the function J^KP/l) sine (a/2)]. 
(D) Behavior of /3 eff //3, where the atomic diffraction patterns are fit with the distribution 
J%(P eS /2). The black solid line is sine (a/2), and the dashed black line corresponds to a 
numerical fit using /3 =4.5. 



The presence of higher diffraction orders leads to population decays that are faster than those 



determined by the two-photon timescale in Eqs. (|9J) and ([TO]). This is shown in inset (C) of Fig. 
[2] for the ±4hk (n = 2) orders, which require the exchange of four photons in the transition 
from Ohk. In general, the presence of many higher momentum orders will complicate analytic 
descriptions of the diffraction dynamics. 

Nevertheless, it is possible to accurately describe the observed dynamics by numerical inte- 
gration of the coupled differential Eqs. ([?]), truncated for unpopulated higher orders beyond n. 
This is shown in Figs. [2] and [3] (A). 

6. Depth calibration of optical lattices 

In experiments with ultracold atoms in optical lattices, the tunneling rate depends exponentially 
on the lattice depth Vq fl9lL for which an accurate knowledge thus is essential. Unlike most 
other methods (cf. 1 19]), applying a "short" Kapitza-Dirac diffraction pulse can conveniently 
reveal the lattice depth in a single-shot measurement. However, the application of well-defined 
pulses that are short enough to be deeply in the Raman-Nath regime, yet strong enough to 
yield significant diffraction, can be technically challenging. The present study quantifies how 
the results of such a calibration measurement need to be corrected if pulses of finite length are 
used. 
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Fig. 3. (A) Breakdown of the analytical description. For each data set, numerical simula- 
tions (using Eqs. ([4])) of the momentum state distribution are shown as solid curves. The 
red crosses are data for the Ohk order, and the open squares (filled dots) are data averaged 
for the ±2hk (±4hk) orders. (B) Comparison of lattice depth calibrations from relatively 
short (8 jUs) and long (50 /is) optical lattice pulses for various intensities of the standing 
light wave. The results for numerically-fit lattice depths agree to within 4% (the solid line 
has a slope of 1). 

In Fig. [3] we present the comparison of experimentally observed momentum state distribu- 
tions to numerical simulations, for three different pulse lengths (8, 25, and 50 /is), as the lattice 
depth is varied from to 65 E r . While the observed momentum state distributions are close to 
the typical Raman-Nath form Eq. ^ for the 8 /is pulses, this is no longer the case for longer 
pulses. For all pulse lengths, numerical simulations agree well, as can be seen in Fig. [3] (A). 
The lattice depth calibrations for the three pulse durations, determined by comparison of exper- 
imental data to numerically simulated patterns, agree to within 4%, as can be seen in Fig.[3](B) 
for the extreme cases of 8 and 50 /is. Despite the complicated form the dynamics take, one 
can thus obtain a reliable lattice depth calibration even with longer pulse durations. We have 
independently verified the method described here with another single- shot calibration method 
for deep lattices (adiabatic lattice rampup, followed by a sudden projection of the ground band 
population onto the ±2hk plane- wave states (24)), for which we found comparable agreement 
to within 10%. 

7. Conclusion 

We have studied Kapitza-Dirac diffraction of a Bose-Einstein condensate in a pulsed standing 
light wave, considering the case that the pulse area remains constant as the pulse duration 
is varied. We find that for sufficiently weak pulses exciting only the ±2fik orders, the usual 
analytical short-pulse prediction continues to hold for longer times, albeit with a modification 
of the apparent modulation depth of the standing wave. We relate this effect to the frequency 
spread of the square-wave pulse, and also draw analogies to the Rabi dynamics of a coupled 
two-state system. Our findings are of practical interest for the calibration of optical lattices in 
ultracold atomic systems, and we show that for a general length and strength pulse, relatively 



simple (neglecting mean-field effects) numerical modeling can be used well outside the Raman- 
Nath regime to accurately determine the lattice depth. 
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